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PROPAGATION OF WAVES OF ACOUSTIC 
FREQUENCIES IN CURVED DUCTS 
by Wojciech Rostafinski 
Lewis Research Center 

SUMMARY 

The propagation of waves of acoustic frequencies in curved ducts is studied for 
the first four modes. The analysis makes use of Bessel functions of the order (n + -|) 
to construct curves of wave number in the duct versus imposed wave number. The 
results apply to ducts of arbitrary width and arbitrary radii of curvature. 

The characteristics of motion in a bend are compared with propagation of waves in 
a straight duct, and important differences in the behavior of waves are noted. 


INTRODUCTION 

Propagation of waves in curved ducts and pipes is characterized by wave patterns 
totally different from those known in straight ducts or in unlimited space (ref. 1). The 
curvilinear boundaries are responsible for the appearance of a continuous standing 
radial wave which in turn affects the transmitted tangential waves. 

However, the method of the basic analysis on long waves (frequency less than 
5 Hz) could not be extrapolated to higher frequencies, and a different analytical approach 
was needed. 

Bessel functions of the order (n + -|) are used to generate characteristic curves of 
the angular wave number versus imposed wave number. By using these curves one 
can obtain the angular wave number for any desired (imposed) frequency and examine 
the motion of waves in bends in detail. 

Analysis of waves in the acoustic frequency range, at the present time up to 
2500 hertz, is of interest in several areas of engineering. Besides problems of prop- 
agation and attenuation of noise in bends it has an application in studies of propagation of 
disturbances through compressor blade rows. 

The analysis which follows concern a stabilized motion in bends with the exclusion 



of the transient effects in junctions between bends and straight lines or any other dis- 
continuities. The physical systems considered consists of a two-dimensional circular 
bend with inner and outer wall radii and Rg, respectively. The walls will be 
assumed perfectly rigid. 


ANALYSIS 

The linearized wave equation in cylindrical coordinates is known to be separable 
in coordinates proper for the boundary, and its general solution may be written in the 
following form: 
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where k = oj/c is the wave number of the imposed motion at the inlet and 

F = J„ (kr) J’ (kR ) - J (kr) J* (kR ) = 
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where the primes indicate differentiation with respect to radius. The function F , 

n 

as shown, was established using the boundary condition d<p/dr = 0 at walls R^ and 

Rg. The integration constants, still undetermined, may be used to match the flow in a 

bend with flow in any ducting upstream and downstream of the bend. In the present 

analysis the term depending on the angular position in the bend and containing these two 

constants will not be evaluated. Only the sustained motion in a bend of infinite length 

(like in a coil) will be examined. Consequently for the far end boundary condition no 

reflection of waves need be considered and we set iC = -D . The term in the bracket 

n n 

reduces to -2D n exp i(-^ n 0). The tangential and the radial vibrational velocities 
v 0 and v . respectively, are as follows: 
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where F' is the derivative of F with respect to r 
u n n 


F’ (kr,kR«) = J' (kr) J (kR.) - J’ (kr) J' (kR.) 
n n n n n 


( 1 ) 


If we set r = Rg = aRj with a = Rg/Rj^ anci ec l ua t e equation (1) to zero, we obtain 

F' (akR., , kR..) = J* (akR.) J ' (kR.) - J* (akR.) J* (kR.) = 0 (2) 

n n n n n 

We are interested in finding the real roots of this characteristic equation, the i^'s. 
They will be fractional numbers, functions of (kR^) and of the parameter a. The 
i^ n 's have the character of an angular wave number and determine the angular phase 
velocity Q - (jJ / v n in curved ducts. Determination of the roots i> n , not available in 
tables, required use of a special procedure. This procedure uses Bessel functions of 
the first and second kinds of the order (n + which are characterized by closed form 
expressions. Using published (ref. 2) expressions for spherical Bessel functions and 
the Raleigh's formula, F^ were calculated for r = aRp = n + \ (n = 0, 1, 2, . . . ) 

and for (kRj), in small increments, ranging from 0. 1 to 20. This study gave a series 
of roots (kRj) satisfying characteristic equation (2). Solution of this inverse problem 
allowed construction of curves and tables which yield, by interpolation, values of 
for any given (kRj). Typical results are shown on figure 1. The graph gives the 
angular wave numbers v's for any arbitrary, imposed, wave number parameter (kRj). 
The calculated curves are for a = R 2 /R 1 = 2.0. The graph indicates that, up to 
(kRj) = 3. 20, only a single mode will be transmitted by a curved duct. At 3. 20 < kR, 

< 6. 30 two modes may be transmitted provided that both modes are present at the 
inlet to the bend. In order to interpolate between v n = (n + ^), the function F' v was 

formed using general expansion of Bessel functions, limited to 15 terms, for the 
arbitrary, noninteger v. The results of calculations are given in table I. To each 
v n correspond two columns of values of (kR^), one calculated from the closed form 
expressions for the Bessel functions, the other using 15 terms of the Bessel series. 

It will be noted that all roots (kR^ of the inverse problem, and consequently, the v 
for the direct problem (when (kR^) are imposed) are calculable by both expressions 
with high accuracy in the basic (zeroth) mode and in the first mode up to = 3. 5. 
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Imposed wave number parameter, kRj 

Figure 1. - Characteristics of motion in a bend for a ■ ■ 2. 


The roots of the second and higher modes are not exactly calculable by the general 
series expression. The effect of the curvature of the duct on the behavior of waves 
will be more apparent by comparing figure 1 to figure 2 , which illustrates the behavior 
of waves in straight ducts. Figure 2 was established using the theory given in several 
texts (refs. 3 and 4). The equation for the velocity potential is .... 
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Imposed wave number, k/7r 


Figure 2. - Characteristics of motion in straight ducts for 2L ■ 0. 2 meter. 


Wave number, v 
5.54.5 3.5 2.5 1.5.5 




Radial velocities 


Figure 3. - Vibrational velocities of zeroth mode in curved duct of 
a * R^Ri = 2 for a range of frequencies. All curves are nor- 
malized to maximum amplitude of 1.0. 


low frequency ( v = 4), the distribution of tangential velocities across the duct’s width 
is very nearly that of a potential vortex, inversely proportional to radius. At a higher 
frequency (v Q = 5-|) the distribution becomes more of the forced vortex type, propor- 
tional to radius. However, the radial velocities exhibit little change in the same range 
of frequencies; a quasi -parabolic distribution remains and only the maximum is 
shifting continuously upwards with frequency. To continue evaluation of the data in 
figures 1 and 2 vibrational velocities have been calculated for four modes and for a 
single angular wave number (v = -|). The results are shown in figure 4. In figure 5 
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Figure 5. - Vibrational velocities in straight duct for one 
frequency and three consecutive modes. 



the same type of data is presented for a straight duct. The vibrational velocities in a 
curved duct exhibit more pronounced changes than in a straight duct. It will be noted 
that the first and the higher modes in a straight duct are characterized by appearance 
of transversal velocities not present in the zeroth mode. In the case of a curved duct, 
the width of a duct or sharpness of the bend have a significant effect on the character- 
istic of the motion of waves. In figure 6, the data of figure 1 are supplemented by 
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Imposed wave number parameter, kRj 

Figure 6. - Characteristics of motion in bends for three bends of different widths. 


curves calculated for a = 1. 5 and 2.5. In figure 7 the data of figure 2 are supplemented 
by data pertaining to a duct only half as wide as the duct of figure 1. Comparing the 
two maps we note that the width of a straight duct has no impact on the character of 
the basic mode. In the curved duct the width of the channel has a significant effect. 

The wave number will change, and the distribution of velocities will be altered. In 
both the straight and the curved ducts, the narrower the duct is the less chances are 
that higher modes will be transmitted. A wide duct admits higher modes much more 
easily. Roots of equatibn-(2) for several values of the wave dumber v n and for two 
values of a are given ifi|able II. They complete the aata of table I. It was found that 
interpolation between dala? of the two tables was satisfactory provided that a curve is 
traced between three points taken from the table sT 
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Imposed wave num^er^w 

Figure 7. - Characteristics of motion in two straight ducts of different width. 
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Finally as an application, let us evaluate motion in a curved duct of R^ = 0. 1 meter 
and Rg = 0. 2 meter (a = 2). We assume that, at the inlet (kR^) = 3.45, frequency 
equals approximately 1900 hertz and that three modes are being generated. From fig- 
ure 1 we learn that the duct will admit only two modes. The angular wave number of 
basic mode will be ^ = 5.5 and that of the first mode is = 1. 74. The distribution 
of the tangential and radial vibrational velocities is shown on figure 8. The vibrational 
velocities of the zeroth mode are positive everywhere. The tangential vibrational 
velocities of the first mode are negative in the outer half of the duct. 




Figure 8. - Vibrational velocities in curved duct of Rj ■ 0. 1 meter with inlet 
(kRj) ■ 3.45. 


CONCLUDING REMARKS 

Propagation of acoustic waves in a frequency range of 5 to 2500 hertz in curved 
ducts was examined. Using an inverse method (calculating the imposed wave number 
parameter (kR^) for an assumed angular wave number v instead of vice versa) cor- 
relation has been obtained between the wave number k of the imposed motion and the 
angular wave numbers in bends v . The results have been tabulated for a range of 
parameters. Analysis indicates that interpolation between the tabulated data by means 
of series expansion of Bessel functions yields accurate results for the zeroth mode 
of motion in the range of arguments (kRj) up to ~5. 0. 
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Analysis showed that distribution of vibrational velocities in curved ducts strongly 
depends on frequency and the angular wave number depends on a = FLj/Rp a measure 
of the bend’s sharpness. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 30, 1972, 

501-24. 
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